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Measurement of spatial resolution of 11 MeV proton radiography* 
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Proton radiography experiment with a Zumbro lens system was carried out on an 11 MeV proton cyclotron. 
The experimental results show that the image blurring is improved markedly. Clear images and good spatial 
resolution of the density step edges are obtained, which is important for hydrotest experiments, and the spatial 


resolution can achieve ~100 um. 
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I. INTRODUCTION 


Proton radiography (PRAD) is an advantageous method 
for hydrotests, with high penetrating power, high detection 
efficiency, small-scattered background, inherent multi-pulse 
capability, and large standoff distances between the object 
and detector [1—4]. The proton-materials interaction renders 
PRAD the capability of resolving higher density objects than 
flash x-ray radiography and identifying the materials. PRAD 
images can be used to evaluate interfaces, the position and 
velocity of the front of shock, and density transformation of 
test samples. The spatial resolution and areal density resolu- 
tion of static objects are benchmark of working status of the 
imaging system. 

In penetrating materials, the protons interact with 
electrons, the nuclear Coulomb field and the nucleus, hence 
the three effects: energy loss, multiple Coulomb scattering 
(MCS) and nuclear scattering. The MCS effect is the main 
cause of the blurring with low energy incident proton beams. 

To eliminate this blemish, Zumbro and Thomas proposed 
a simple, elegant magnetic lens system to focus the scattered 
protons. Then, magnetic lens systems were designed for pro- 
ton radiography experiments on the 800 MeV accelerator at 
Los Alamos Neutron Scattering Center (LANSCE) [5] and 
on the Alternating Gradient Synchrotron (AGS, 24 GeV pro- 
ton beams) at Brookhaven National Laboratory [6]. The ex- 
perimental results proved that this kind of lens system could 
reduce the image blurring remarkably. In their experiments, 
the spatial resolution was about 100 um. 

Protons in tens of MeV can penetrate samples in thickness- 
es of hundreds of microns, making the proton beams become 
a proper tool for imaging sub-mm samples. Also, PRAD may 
be used for imaging medical and biological samples. How- 
ever, for low energy PRAD, the spatial resolution is a key 
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parameter to get the sample information. In this paper, we re- 
port the proton radiography with unit magnification Zumbro 
magnetic lens system on an 11 MeV cyclotron. By analyzing 
the radiographed images of static objects, the spatial resolu- 
tion of magnetic lens system is about 100 um, which is an en- 
couraging result for developing the technique for hydrotesting 
experiment. 


I. METHODS 
A. The chromatic cancellation of Zumbro lens system 


The Zumbro lens system is a point-to-point image system 
with the unit magnification [7, 8]. Along the x and y direc- 
tion, the transport transfer matrices are —I’s. The system is 
composed of two identical doublets of quadrupole, which are 
reflection symmetric. 


B. Blurring analysis 


In radiography, the Zumbro lens system, the radiographed 
objects, the scintillator and the optical recording system all 
contribute to the image blurring. Blurring contributions from 
each part are considered as follows. 

For the Zumbro lens system, the position-dependent chro- 
matic aberrations, and the aberrations and blurring brought 
by the transport of proton beams, can be eliminated because 
that the lens works on the point-to-point radiography. Howev- 
er, for different spots of the sample, the chromatic aberration 
cannot be eliminated. The chromatic aberration blurring of 
different positions of the sample is calculated as [7] 


Biens = T126906, (1) 


where ĝo is deviation angle of the MCS, 6 = Ap/po is the 
momentum deviation, and T}2¢ is element of transport matrix 
in TRANSPORT notation, which represents the second order 
chromatic aberration. 
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Fig. 1. (Color online) Schematics of the radiography beamline. 


Before exiting the object, the protons have undergone mil- 
lions of Coulomb scatterings, and angle distribution of the 
penetrating protons can be treated approximately as Gaussian 
distribution. The root-mean-square (rms) width, in milliradi- 
ans, of MCS exit cone is 


Oo = [13.6/(p8c)](L/Lr)'/*{1 + [log(L/Lr)]/9}, (2) 


where p is the beam momentum in GeV/c, (c is the ve- 
locity of the proton, L is the thickness of the object, and 
Lg x Z? /A is the radiation length of the material [6]. So 
the blurring coming from the object can be considered as [9] 


Boy = 37"? L00. (3) 


The blurring from the detection system, i.e. the scintillator 
and the downstream optical image system (the optical lens 
and the CCD camera) is denoted as Baet. Therefore, the total 
blurring can be calculated as, 


1/2 
Brot FWHM = 2.36 (Boj + Bras + Bie) f ; (4) 


where the coefficient 2.36 is obtained from the Gaussian dis- 
tribution. This total blurring is the theoretical FWHM (Full 
Width at Half Maximum) of the edge differential image. 


C. Experimental setup 


The PRAD experiment for spatial resolution of the Zumbro 
lens system with static objects was carried out on the 11 MeV 
proton compact cyclotron [10], which provides average beam 
current of 50 uA. The collimated proton beam on the scintil- 
lator detectors was integrated from thousands of micro pulses 
in ten seconds or longer. 

Figure 1 is a layout of the low energy PRAD imaging 
beamline [11]. It consists of a matching section, the mag- 
netic lens system and the recording system. The matching 
section includes two correction magnets and two quadrupoles 
to make the beam matching to the magnetic lens system. The 
recording system includes a LSO scintillator at the image 
plane and a high resolution CCD camera about 15 cm down- 
stream of the LSO scintillator. 

The object was assembled by aluminium foils of differ- 
ent thicknesses and charcoal layer coated by a laser printer. 
Thickness of the charcoal layer was determined as follows: A 
50mm x 50mm x 6 um Al foil of known weight was coated 


with a charcoal of known density. It was weighed by an elec- 
tric balance of 10 ug precision. Then, the averaged thickness 
of charcoal was determined as 4 um. 

A ¢6mm collimator was fixed on the Fourier plane, sort- 
ing the angle of scattered protons. Its maximum picking an- 
gle is 4.06 mrad in our experiment. The field of view (FOV) 
is larger than 630 mm. The scintillator is a 638mm LSO 
(Lu, SiO,) screen, and the camera is a high quality industrial 
CCD camera. 


HI. EXPERIMENTAL DATA 


Parameters of the imaging system are known or calculated. 
For example, Ti26 = 2929mm is determined by the design 
of the Zumbro magnetic lens system, and the total blurring of 
the image beamline with an object of Al foil in thickness of 
16 um is estimated as follows: Because the object is too thin 
to affect the momentum deviation, the 6,m; should be momen- 
tum deviation of the cyclotron, which is about 0.003. A half 
of momentum deviation comes from the cyclotron and the 
other half comes from the diffuser. From Eq. (2), the MCS 
angle is 99 ~ 5.6 mrad without the collimators. Therefore, 
the blurring from lens system is Biens = 82 um. With the col- 
limator, the rms of MCS angle changes. After the first dou- 
blet, the exiting protons are collected as a spot on the Fourier 
plane, the collimator on the spot position cuts off the pro- 
tons with a smaller diameter, the maxima angle changes to 
the collimated angle. And on the scintillator, the angle @ is 
3.38 mrad, and the blurring is Biens = 29.73 um. For an Al 
foil of 9 um thick, the blurring with the collimator can be es- 
timated at Biens = 29.70 um. 

According to Eq. (3), the amount of 0.064 um and 0.04 um 
body blurring (,pj;) correspond to the 16 um thick Al foil and 
the 9 um thick Al foil, respectively. 

From SRIM calculations, for 11 MeV protons travelling a 
0.5mm thick LSO scintillator, the lateral straggling is about 
36 um. The blurring Bg, is about 36.9 um. So in our experi- 
ment, from Eq. (4), the total blurring (FWHM) for Al foil of 
16 and 9 um in thickness is Biotpwom © 112 um. 

Figure 2 is the transmission image and edge position image 
of the stepped sample with a ¢6 mm collimator on the Fourier 
plane. The density step positions in Fig. 2(b) are identified as 
those points where the gradient of image is the maximum. To 
get more precise position of the edge, we analyzed the red re- 
gion in Fig. 2(a). By averaging the region along Y direction, 
and differentiating it along X direction, we had Fig. 2(c), and 
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Fig. 2. (Color online) Proton transmission images of Al foils in 
thickness of 0, 6, 9, 15, 16, 22, 24 and 31 um (a), the edge posi- 
tions processed by MATLAB (b), and Gaussian fitting of differential 
vertical gray scale along the Fig. 2(a) red-line region (c). 
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Fig. 3. (Color online) The transmission image of reticle radiography. 


then made a Gaussian fitting with the differential peaks. The 
maxima of the differential peaks and the FWHM values were 
obtained. Comparing with the edge image, the peak is at the 
same position. FWHM of the peak is about 100 um, which 
means that the spatial resolution of total imaging system is 
around 100 um. The distance between two edges in Fig. 2(b) 
is 5 320 mm, while in Fig. 2(c), itis 5 308 mm. For the simple 
configuration, searching for the interface position of the dif- 
ferential using MatLab toolbox can provide cursory positions 
of the interface edges. 

To obtain more direct information of the spatial resolution, 
check the spatial resolving power, and confirm the feasibility 
of searching for the interface position, reticle samples were 
radiographed. Fig. 3 shows transmission images of reticles of 
0.5, 1.0, 2.5 and 5.01p/mm prepared with a laser printer on 
a 50mm x 50mm x 6um Al foil. We analyzed the intensi- 
ty distribution and the derivative of the intensity distribution 
along the rectangle region. 

In Fig. 3, the 2.51p/mm reticle can be distinguished clear- 
ly, and the 5lp/mm reticles can also be distinguished by 
eyes, so 100 um spatial resolution of the whole system can 
be achieved. 

Comparing the results above, there is an inconsistency be- 
tween the experimental data and the theoretical image blur- 
ring. They may come from the complex environment such as 
the vibration of vacuum pumps, and the jitter of proton beam. 
To obtain clearer radiography, efforts will be made for further 
understanding of the principle and working status of Zumbro 
lens system. 


IV. CONCLUSION 


We have performed the spatial resolution experiment of 
proton radiography on the 11 MeV cyclotron with a unit mag- 
nification Zumbro lens system. By analyzing the transmission 
image of reticle and edges image of Al foils, the data show 
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that the imaging system can achieve the resolving power of 
100 um. 

The results of this experiment are encouraging. They prove 
that the static radiography on the low current cyclotron is fea- 
sible. However, it is still a lot of work to improve spatial 
resolution of the system. For example, under bombardmen- 
t of 11 MeV protons hit on, the scintillation occurs not in a 
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